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Abstract. A framework for statistical-mechanical analysis of quantum Hamiltonians 
' is introduced. The approach is based upon a gradient flow equation in the space of 

, Hamiltonians such that the eigenvectors of the initial Hamiltonian evolve toward those 

ly^ ' of the reference Hamiltonian. The nonlinear double-bracket equation governing the 

flow is such that the eigenvalues of the initial Hamiltonian remain unperturbed. The 
space of Hamiltonians is foliated by compact invariant subspaces, which permits the 
• construction of statistical distributions over the Hamiltonians. In two dimensions, an 

explicit dynamical model is introduced, wherein the density function on the space of 
^ ^ Hamiltonians approaches an equilibrium state characterised by the canonical ensemble. 

■ This is used to compute quenched and annealed averages of quantum observables. 
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^ , In the conventional approach to statistical mechanics the Hamiltonian of the system 

^ ! under consideration is held fixed. If the system is in equilibrium with a heat bath, 
00 
O . 

QO , interactions with the bath. The equilibrium distribution over the state space of the 
system (configuration space of a classical spin system, classical phase space, or the space 
of pure quantum states) is then established. However, in some cases — as in amorphous 
alloys — the Hamiltonian need not be fixed, and may even fiuctuate owing to thermal or 
5^ \ other intrinsic sources. Observable effects arising from such Hamiltonians may even be 
significant in the quantum domain. 

The purpose of the present paper is to introduce a theoretical framework for an 
equilibrium theory of Hamiltonians. The fact that parameters or matrix elements of 
the Hamiltonian themselves are subject to random fiuctuations for some systems has 
long been recognised in the literature of spin glass [1] or random matrix theory [2]. The 
novel idea introduced here, as distinguished from that considered in the theory of spin 
glass or random matrices, is the construction of equilibrium distributions over invariant 
subspaces of the space of quantum Hamiltonians by using a gradient fiow equation on 
the space of Hermitian matrices. 

In classical statistical mechanics the notion of a gradient fiow plays an important 
role in describing the approach to equilibrium: A system immersed in a heat bath 
naturally tends to release its energy into the environment and thus approach its 
minimum energy state, and this tendency is characterised by a Hamiltonian gradient 
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flow. An equilibrium state is attained when this flow is on the average counterbalanced 
by thermal noise, where the magnitude of the noise is determined by the temperature 
of the bath. Accordingly, we shall introduce a gradient flow equation on the space of 
Hamiltonians such that the eigenstates of an arbitrary initial Hamiltonian Hq at time 
t = tend toward alignment with those of a reference Hamiltonian, denoted by G. 
Thus, G plays the role of the 'fixed' Hamiltonian in conventional quantum statistical 
mechanics. The eigenstates of Ht thus evolve toward those of G under the flow. By 
introducing of a suitable noise term, we then characterise the approach to an equilibrium 
distribution. 

The paper is organised as follows. The key results concerning the properties of the 
double-bracket equation that generates the gradient flow are summarised first in the 
Proposition. The notion of a double-bracket flow was first introduced by Landau and 
Lifshitz in the context of characterising dispersions in magnetism [3]. In its 'modern 
form' it was introduced by Brockett P] and has been successfully applied to many 
areas, such as optimal control, linear programming, sorting algorithms, and dissipative 
systems. Although some assertions of the Proposition are valid in all dimensions, we 
shall analyse only the two-dimensional case in full detail. We subsequently construct 
an explicit model for the 'equilibrisation' of 2 x 2 quantum Hamiltonians, such that the 
stationary state is given by the canonical distribution. The resulting statistical theory 
of quantum Hamiltonians can be extended to a modification of quantum statistical 
mechanics. In particular, we work out the quenched and annealed averages of quantum 
observables. We conclude by indicating how the analysis can be extended to higher 
dimensions. 

Proposition. Let Ht and G be arbitrary 2x2 Hermitian matrices, where Ht is 
time dependent and G is fixed. Let Ht satisfy the double-bracket evolution equation 



with initial condition Hq. Then the evolution ([T]) is isospectral, i.e. the eigenvalues of Hq 
are preserved under ([1]), and \im.t^oo[HtiG] = 0. Furthermore, the space of Hermitian 
Hamiltonians is foliated by a family of invariant 2- spheres L, and ([1]) induces a gradient 
flow on each C 

We remark that in terms of the Hermitian operator X = i [H, G] the double- 
bracket evolution ([T]) can be rewritten as dH = iX[H, X]dt, which formally is just 
the Heisenberg equation of motion. However, owing to the if-dependence of X the 
evolution is nonunitary. We also note that in units h = 1 the parameter A has dimension 
[Energy]"^. The Hamiltonians Hq and G are both assumed nondegenerate; otherwise, 
if at least one of the Hamiltonians is degenerate, then Hq is a fixed point of the flow. 
We now proceed to establish the Proposition. 

The fact that equation ([1]) asymptotically drives Ht toward [Ht, G] = 0, irrespective 
of the dimensionality of the matrices, follows from the relation 



dHt 





(1) 



dt 



TjHHt-Gf 



2tT{[G,Ht]^[G,Ht]) <0 



(2) 
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where the equahty is attained if and only if [Ht, G] = 0. To see that ([T]) defines an 
isospectral flow (which is also valid irrespective of the dimensionality of the matrices) 
we note that the right side of ([I]) can be written in the form Ad(e~'*"^ifte^'^^)/ds|s=o. 
The isospectral property then follows from the relation det(e~'*^ifie'*"^ — El) = 
det{Ht — El). To prove that the orbit of the flow for a given initial value Hq lies on a 
two-sphere C (which is isomorphic to the space of pure states for a two- level system), 
and that ([1]) defines a gradient flow on £, we shall solve ([1]) explicitly for the case of 
2x2 Hermitian matrices. 

Let the 2x2 Hamiltonian be represented in terms of the Pauli matrices as 

Ht = ^Utl + i z/cr-rit, (3) 
where rij = (xj,yt,Zf). Similarly for the reference Hamiltonian G we write 

G = lvl + liJ(T-g (4) 
for a unit vector g. Bearing in mind the relations 

ii oc tT[Ht,X] = and [Ht,G] = |ii//icr-(nt x g) (5) 
we find that ([T]) reduces to 



dnt 



dt 

where u = Xufi. From 
d(nf ■ nt 



untx {nt X g), (6) 



oc ■ {nt X {nt x g)) = (7) 

we see that the norm of remains constant under ([6]). Without loss of generality we 
work with the basis in which G is diagonal, and choose g = (0, 0, 1). In terms of the usual 
spherical parametrisation in the G-basis we have = (sin 0^ cos 0^, sin 6'^ sin 0^, cos 6'*). 
Therefore, ([6]) reduces to: 

dt = ujsinOt and <pt = 0. (8) 
Solving these, we obtain 

cos 6*4 = tanh(co — ut) and (pt = (po, (9) 
where Cq = tanh~^(cos6'o) and 0o are initial values. The solution Ht to ([T]) is thus 

i( Uo — i'tanh{ujt — Co) sech(ci;t — Cq) e~'*' \ 
^ I iysech.{ujt — Cq) e"^" 'Uo + i^tanh(co't — Cq) / 

A straightforward calculation shows that the eigenvalues of Ht are time-independent, 
and that 

limi7, = |f ° V (11) 

Thus, the Hamiltonian is asymptotically diagonalised in the G-basis. Observe that tiHt 
and detHt are conserved quantities. Therefore, the flow induced by ([6]) for fixed initial 
values Mo and |no| is confined to a two-sphere £, which can be identified with the state 




Figure 1. Flow on the sphere C. The vector field generated by tlic unitarily-modified 
gradient-flow is plotted. The first term in generates a rotation around the G- 
axis, while the second term generates geodesic flows toward the south pole. The axis 
no of the initial Hamiltonian Hq spirals around the G-axis g and is asymptotically 
aligned with the latter. 



space of a two- level system (i.e. the complex projective line). Since uq and |no| are 
constant, in what follows we shall fix these two variables and focus our attention upon 
the associated sphere £ parameterised by the dynamical coordinates {6t,4't)- 
The fact that (E]) defines a gradient flow 

dx'^ = -\ Az/^"''VbG(x)dt (12) 

on C, where we use local coordinates {x^,x^) = {9,(f)) on the sphere, can be seen as 
follows. First, in terms of these coordinates the inverse metric on the sphere is 



9^' = —2-n\ : 1- (13) 




We define a function G{x) on the sphere £ as follows: 

G{9,^) = l{v + ficos9). (14) 

This is obtained by taking the 'expectation' of reference Hamiltonian G in a pure state 
corresponding to the point (6*, 0) on C Then, a short calculation using (IT^ . (IT^ . and 
([T^ shows that the dynamical equations ([8]) correspond to the gradient flow ([T2|) . 

We remark, incidentally, that the dynamical equation ([T]) can be modified to include 
a unitary term: 

^ = -i[H„G]-X[H„[H„G]], (15) 

without greatly affecting its physical characteristics. In the 2x2 example considered 
here, the only change occurs in the phase, so that instead of (pt = 4>o we have 0t = 4)0+ fit. 
Thus, according to ([T]) the eigenvectors of Hf evolve 'straight' toward those of G (i.e.. 
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along geodesies), whereas under (IT^ they 'spiral' toward those of G. This is illustrated 
in Figure [1] where we plot the vector field on the sphere defined by the dynamical 
equation ( ITSi) . 

We also note that the diagonalisation property of the double-bracket evolution 
dl]) has been applied to the analysis of Toda lattices [5], dispersions in the Euler- 
Poincare equations [6], couplings in photorefractive media [7], and fiow equations in 
renormalisation group [8]. Although here we consider the case in which G is fixed, it 
is possible to vary G in time (that is, vary the direction of g). Then the dynamical 
equation can be used to characterise quantum control (cf. [9J for a related idea). 
In this context it would be interesting to investigate the role of geometric phases for 
observables, when the control Hamiltonian G is varied along a loop in C. 

Having defined a natural gradient fiow on the invariant 2-spheres foliating the space 
of Hamiltonians, we now consider a dynamical model on a given sphere C such that an 
arbitrary initial Hamiltonian Ho evolves — according to ([1]) — towards the reference frame 
determined by G, but at the same time is randomly perturbed in all directions in C by 
a pair of independent Brownsian motions. The dynamical model, in particular, will 
possess the following properties: (i) the eigenvalues of Ht remain constant in time, 
and (ii) the probability distribution over C evolves toward an equilibrium distribution 
characterised by the standard canonical density function. Although rather elaborate, 
this model can be treated analytically by identifying any given surface of the foliation 
with the space of pure states of a two-level system, which permits application of the 
model for the thermalisation of quantum states introduced in |10] . 

We consider first a stochastic differential equation of the form 

dx"" = /iMt + na^dW^ (16) 

on C (viewed as a real two-sphere). Here k is a constant, the drift /i" is a vector field 
on C, and the vectors {o",?}j=i 2 constitute an orthonormal basis in the tangent space 
of £ such that g""^ = afa'jS'^^ and afajgab = Sij. We note that dx" is the covariant 
Ito differential [11], and that the standard 2-dimensional Wiener process {W^} satisfies 
dW^dW^ = S^Mt. By straightforward calculation one verifies [10] that the density 
function pt(x) on £ associated with the stochastic evolution (ITGll satisfies the Fokker- 
Planck equation 

^pt(x) = -V.(/i>t) + (17) 

For our model we require that the drift vector /x" represent the double-bracket gradient 
fiow ([T]). This is achieved by choosing p" = — | k^AV^G, where = z/. Then it follows 
from a theorem of Zeeman [12] that there exists a unique stationary solution to f|T7|) . 
given by the canonical density 

jj,exp[—XG{x))dV 

To illustrate these results in more explicit terms we consider a system consisting 
of a single spin-i particle immersed in an external magnetic field. The Hamiltonian is 
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then H = — B ■ S, where B denotes the field and S the spin vector. The direction of 
the field B, however, is subject to fiuctuations around its stable direction, specified by 
G (directed along the z-axis). Calculating the orthonormal basis on the sphere, we 
obtain the stochastic equations for the variables {0,(j)): 

dOt = u sin ^tdt + v^(dW; + dW^) 

d0, = --^v^(dW^,i-diy,2). ^ ^> 



The associated Fokker-Planck equation reads 

p = -uj{cose + sinede)p + 2u{dj+ -tr-sdl)p, (20) 

where dg = d/dO and = d/dcp. The asymptotic solution is the following canonical 
density function: 

Direct substitution shows that ( 12T1) is the stationary solution to ( l20l) . It follows from 
( !2T|) and the use of the spherical (Fubini-Study) volume element dV = j sin 6'd6'd0 that 
the equilibrium mean Hamiltonian is 

where {cos9)x = 2/A/i— 1/ tanh(|A/i). We may regard the parameter A as representing 
the 'inverse temperature' for the Hamiltonian: if the noise level is high (A <^ 1), then 
the direction of the external field B on the average lies close to the xy-p\a.ne so that 
(cos 6*) A — 0; whereas if the noise level is low A ^ 1, then the field B on the average is 
parallel to the 2;-axis and we have {cos6)\ ~ —1. 

Let us now consider how the statistical theory of Hamiltonians presented here can 
be applied to quantum statistical mechanics. In this context it is natural to borrow 
ideas from the spin glass literature [Tj. We may take the averaged Hamiltonian {H)x as 
the starting point of the analysis — this gives the analogue of an annealed average. In 
this regime the expectation of an observable O is given by 

Such an expectation, however, will involve the use of the averaged Hamiltonian {H)x 
whose eigenvalues differ from those of H. Alternatively, we may use the 'unaveraged' 
Hamiltonian to compute the thermal expectation of an observable O, regarded as a 
function on a specified invariant surface in the above-described foliation of the space of 
Hamiltonians, and then take its average — this gives the analogue of a quenched average: 

A short calculation shows that the canonical quenched average of the Hamiltonian G is 
(G)g = i/.tanh(l/3.)(^^^-A^, (25) 
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whereas the canonical annealed average of G is 



{G)a = ^/itanh 



2 



(26) 



These averages are plotted in Figure [2l These results suggest a new line of studies on 
the extended quantum statistical mechanics of disordered systems. 

The explicit analysis presented above is for the most part confined to 2 x 2 
systems. In higher dimensions, the double-bracket evolution equation ([T]) still defines 
an isospectral gradient flow in the space of Hamiltonians. Thus, the procedure for a 
statistical analysis of Hamiltonians as outlined above is naturally extendable to higher 
dimensions. However, in higher dimensions the equivalence of the Schrodinger and 
Heisenberg picture for the nonunitary motion breaks down (that is to say, the 
generic surface foliating the space of Hamiltonians is not isomorphic to the associated 
space of pure states). Instead, in higher dimensions, the relevant foliation consists 
of certain subspaces of higher-dimensional spheres. Nevertheless, there exist unitary- 
invariant measures on these spaces, which can be used to formulate the theory in an 
analogous manner. In particular, the equilibrium state resulting from the thermalisation 
dynamics remains canonical in the sense that it is proportional to the canonical density 
exp(— Atr(Gif)) just as in the 2x2 example (cf. [I3]). The remaining open problem is 
the precise geometrical description of the relevant gradient flows in higher dimensions, 
and the specification of the associated measures to calculate partition functions. 
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